A Tale of Two Electrons: Correlation at High Density 



Pierre-Frangois LootQ and Peter M. W. GilQ 
Research School of Chemistry, Australian National University, 
Canberra, Australian Capital Territory, 0200, Australia 
(Dated: August 31, 2010) 

We review our recent progress in the determination of the high-density correlation energy E c in 
two-electron systems. Several two-electron systems are considered, such as the well known helium- 
like ions (helium), and the Hooke's law atom (hookium). We also present results regarding two 
electrons on the surface of a sphere (spherium), and two electrons trapped in a spherical box (bal- 
lium). We also show that, in the large-dimension limit, the high-density correlation energy of two 
opposite-spin electrons interacting via a Coulomb potential is given by E c ~ — l/(8D 2 ) for any 
radial external potential V(r), where D is the dimensionality of the space. This result explains the 
similarity of E c in the previous two-electron systems for D — 3. 
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I. INTRODUCTION 

The Hartree-Fock (HF) approximation ignores the cor- 
relation between electrons, but gives roughly 99% of the 
total electronic energy [T]. Moreover, it is often accurate 
for the prediction of molecular structure 0, computa- 
tionally cheap and can be applied to large systems, es- 
pecially within local (linear-scaling) strategies [5HTU] . To 
reduce the computational cost still further, various nu- 
merical techniques have been developed including, for ex- 
ample, density fitting (or resolution of the identity) |1H - 
ITfj] , pseudospectral and Cholesky decomposition [T?H2"5] , 
dual basis methods [2~4"H2"9"] . and both attenuation [5tJll3"T] 
and resolution [32 -35 of the Coulomb operator. 

Unfortunately, the part of the energy which the HF 
approximation ignores can have important chemical ef- 
fects and this is particularly true when bonds are formed 
and/or broken. Consequently, realistic model chemistries 
require a satisfactory treatment of electronic correlation. 

The concept of electron correlation was introduced by 
Wigner and defined as 



E C = E — E] 



HF 



(1) 



by Lowdin [37], where E is the exact non-relativistic en- 
ergy. Feynman refered to E c as the "stupidity energy" 
[58] because of the difficulties associated with its charac- 
terization in large systems. 

Even though it is a formidable challenge to determine 
the correlation energy accurately, even in simple systems, 
recent heroic calculations on the helium atom [59TH2"] 
have demonstrated how near-exact results can be found. 
Indeed, this elementary chemical system has been com- 
pared to the number it by Charles Schwartz [43:: "For 
thousands of years mathematicians have enjoyed compet- 
ing with one other to compute ever more digits of the 
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number n. Among modern physicists, a close analogy 
is computation of the ground state energy of the helium 
atom, begun 75 years ago by E. A. Hylleraas." 

Although E c in the helium atom is now known very ac- 
curately, certain correlation effects remain incompletely 
understood and, for example, even the Coulomb hole [U] 
itself is more subtle than one might imagine. The pri- 
mary effect of correlation is to decrease the likelihood of 
finding the two electrons close together and increase the 
probability of their being far apart. However, accurate 
calculations have revealed the existence of a secondary 
Coulomb hole, implying that correlation also brings dis- 
tant electrons closer together [45] . The same observation 
has been made in the H 2 molecule by Per et al. [H] and 
it appears that secondary (or long-range) Coulomb holes 
may be ubiquitous in two-electron systems [47] . 

In order to get benchmark results for the development 
of intracule functional theory (IFT) 48-53], we have re- 
cently initiated an exhaustive study of two-electron sys- 
tems [55]. In the present Frontier Article, we review our 
recent progress in the determination of the correlation 
energy in various high-density two-electron systems: the 
helium-like ions (Sec. II A), two electrons on the surface 
of a sphere (Sec. |HB ), the Hooke's law atom (Sec. 
and two electrons trapped in a spherical box (Sec. 



IIC 



II D 



It is reasonable to ask whether an understanding of 
the high-density regime is relevant to normal chemical 
systems but it turns out that most of the high-density 
behaviour of electrons is surprisingly similar to that at 
typical atomic and molecular electron densities. Much 
can be learned about the languid waltz of a pair of elec- 
trons in a covalent bond from their frenetic jig in the 
high-density limit. Moreover, it has led to an under- 
standing of key systems, such as the uniform electron 
gas [SHI 153 1 which form the cornerstone of the popular 
local density approximation in solid-state physics [55] , 
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We also show (Sec. Ill) that, in the large-dimension 
limit, the high-density correlation energy of two elec- 
trons is given by a simple universal rule which is in- 



2 




FIG. 1. The basic coordinates of a two-electron system. 



dependent of the external confining potential. Just as 
one learns about interacting systems by studying non- 
interacting ones and then introducing the interaction per- 
turbatively, one can understand our three-dimensional 
world by studying high-dimensional analogues and intro- 
ducing dimension-reduction perturbatively. 

In this study, we confine our attention to the S ground 
states of two-electron systems. This allows us to ignore 
the spin coordinates and focus on the spatial part of the 
wave function. Atomic units are used throughout. 



II. 



HIGH-DENSITY LIMIT 



For two electrons confined in a spherically-symmetric 
external potential V(r), the Hamiltonian is 



H 



Vf 
2 



V 2 

v 2 



+ V(n) + V{r 2 ) + 



1 

T\l' 



(2) 



where the first two terms represent the kinetic energy of 
the electrons and l/ri2 = 1/ \r\ — r%\ is the Coulomb op- 
erator (Fig. [I]) . After a suitable scaling of the coordinates 
and energy |59 , 60 , the Hamiltonian can be recast as 



H = 



V 2 V 2 11 

2 2 Z ri 2 



where Z measures the confinement strength. Equation 
([3]) is well poised for a perturbation treatment in which 
the zcroth- and first-order Hamiltonians are 



(0) 

2 > 



1 



ri2 



and the one-electron Hamiltonian is given by 



1(0) 



1 



V? + V(ri). 



(4) 



(5) 



The zeroth-order wave function satisfies the eigenequa- 
tion 



ff (0) *o (r 1) r 2 )=^°^o (n,r 2 ). 
and the zeroth- and first-order energies are 



(6) 

(7) 
(8) 



Following Hylleraas [6T] , we can use perturbation theory 
to expand both the exact [61] and Hartree-Fock (HF) [62] 
energies as series in 1/Z, yielding 

E(Z,D,V) = E {0) (D,V)Z 2 + E {1) (D,V)Z 

+ E (2) (D,V) + E (3 \D,V)Z- 1 + ..., (9) 



and 



E HF (Z, D, V) = E { ° F (D, V)Z l + E^(D, V)Z 

+ E {2 p (D, V) + Egl(D, V)Z- 1 + ..., (10) 

where D is the dimensionality of the space. It is straight- 
forward to show that 

E<®(D,V) = E™(D,V), 



EW(D,V)=E§>(D,V) 



(11) 

-0mm ( 12 ) 
and therefore, in the high-density (large-Z) limit, we find 
E^(D,V) = lim E C (Z,D,V) 



lim [E(Z,D,V)-E UF (Z,D,V)} 



Z->c 



E&>(D,V)-E™(D,V) 



(13) 



Helium 



As a first example, we consider the D-dimensional 
helium-like ions (He) [51] [53] where the electrons move 
in the Coulomb field of a nucleus with charge Z, i.e. 



V(r) = -- 



From the foregoing Section, we have 

r d 2 D - 1 d 

1 

dr 2 



'0 



r dr 

and the zeroth-order wave function is 

4 D r 2(n+r a ) 



*o(n,r 2 ) 



(D- 1) D T(D) 



exp 



D - 1 



The E^ and E^ values are given by 



£<°>(Z>,He) = - 
E^{D,Uc) 



(D- l) 2 ' 
4 r(D + $T(Bp) 



(D-1) 2 r(£» + i)r(f; 



(14) 
(15) 
(16) 

(17) 
(18) 



where T is the Gamma function [M] . 

To compute the second-order energy E^ 2 \ we use the 
Hylleraas method |61[ . adopting the length and energy 
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TABLE I. Second-order energies and limiting correlation energies in two-electron systems. 



System 


D = 2 


D = 3 


D = 4 


D = 5 


D = 6 


D = 7 






Second-order exact energies, 


-E (2) (D,V), from (9) 




Helium 


0.632740 


0.157666 


0.070044 


0.039395 


0.025208 


0.017501 


Spherium 


0.227411 


0.047637 


0.019181 


0.010139 


0.006220 


0.004189 


Hookium 


0.345655 


0.077891 


0.032763 


0.017821 


0.011153 


0.007622 


Ballium 


0.057959 


0.014442 


0.006194 


0.003333 


0.002037 


0.001352 






Second-order HF energies, - 


■E™(D,V), from ( 


Tof 




Helium 


0.412607 


0.111003 


0.051111 


0.029338 


0.019020 


0.013325 


Spherium 




















Hookium 


0.106014 


0.028188 


0.012904 


0.007382 


0.004776 


0.003342 


Ballium 


0.324120 


0.069618 


0.028107 


0.014770 


0.008977 


0.005983 






Limiting correlation energies 


-E {2) (D,V), from 






Helium 


0.220133 


0.046663 


0.018933 


0.010057 


0.006188 


0.004176 


Spherium 


0.227411 


0.047637 


0.019181 


0.010139 


0.006220 


0.004189 


Hookium 


0.239641 


0.049703 


0.019860 


0.010439 


0.006376 


0.004280 


Ballium 


0.266161 


0.055176 


0.021913 


0.011437 


0.006940 


0.004631 



scaling of Herrick and Stillinger [65] and employing the 
conventional Hylleraas basis functions [5T] 

1> u = Tpnlm = s n t 2l u m exp(-s/2), (19) 

where uj — (n, I, m) are non-negative integers and 

s = n + r 2 , t = r\-ri, u = r 12 , (20) 

are the conventional Hylleraas coordinates. The second- 
order energy, which minimizes the Hylleraas functional, 
is then given by 

£: (2) (L>,Hc) = ~b T A ^b, 



where 
A, 



D - 1 



(21) 

(22) 
(23) 



In (22) and (23), T, L, S and U are the kinetic, electron- 



nucleus, overlap and repulsion matrices, respectively, and 
are defined by 

+ s(u 2 - t 2 ) (d a ip Ul d u ip U2 + duip^dsip^) 
+t(s 2 - t 2 ) (9tVa>i9«V , a>a + du^dtip^)] dr 



and 



s„ 



IpunlpuiadT, 



u. 



T\2 



(24) 

(25) 
(26) 
(27) 



with the volume element and domain of integration [65j 
dr = u(s 2 -t 2 ) J^dsdtdu, (28) 
J = {s 2 u 2 ) (u 2 t 2 ) , (29) 



dr = ds du dt. (30) 
Jo Jo Jo 

All the matrix elements can be obtained in closed form 
using the general formula 



r(n + I + m + 2D - 3) 
,'m + l D-l 

= B 



where 



B 



B(m, n) = 



s n t l u m e - s J^dsdtdu 



l + m + D-1 D-l 



r(m)r(n) 



(31) 



(32) 



r(m + n) 

is the beta function [M] . The E 2 value for D — 3 has been 
studied in great detail [551 [57J , but the only other helium 
whose E 2 value has been reported is 5-helium |65j and 
this value was obtained by exploiting interdimensional 
degeneracies [55] , Numerical values of for 2 < D < 7 
are given in Table |T) 

(2) 

£jjp values can be found by generalizing the Byers- 
Brown-Hirschfelder equations [59] to obtain 

r w (r? dr 

J r D -i* (r,r) ' 
W(r) = 2 / [J{x) - V (x, x) x°- x dx 



E 



(2) 
HF 



J(r) 



max(r, x) 



3-D 1 D 



2' 2 



,a 



(33) 

(34) 

D - l dx. 
(35) 
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where a 



i(x,r) 



and F is the Gauss hypergeometric 



max(a:,r) 

function [M] , 

E^ F (D,He) has been reported for D = 3 by Linder- 
berg [62] and Eq. ( 33 ) yields expressions such as 



E 



(2), 



9 



HF (3,He) = — ln- 



13 

432' 



£7<£(5,He) 
4 2 F(7,He) 



903 , 3 

in - 

1024 4 



35 213 



5 643101 j 3 
' 204 800 n 4 



124416 

640 149 405 049 



(36) 
(37) 
(38) 



80 621568 000 
Numerical values for 2 < D < 7 are shown in Table |TJ 

B. Spherium 

Spherium (Sp) consists of two electrons, interacting 
via a Coulomb potential but constrained to remain on 
the surface of a sphere of radius R = 1/Z [T0TIT3"] . This 
model was introduced by Berry and co-workers [74Tf77] 
who used it to understand both weakly and strongly cor- 
related systems, such as the ground and excited states 
of the helium atom, and also to suggest an alternating 
version of Hund's rule [78] . Seidl studied this system 
in the context of density functional theory [TH] to test 
the ISI (interaction-strength interpolation) model [50] . 
More recently, we have performed a comprehensive study 
of the spherium ground state, using electronic structure 
methods ranging from HF theory to explicitly correlated 
treatments [70] . 

In this Section, we consider Z3-spherium, the gener- 
alization in which the two electrons are trapped on a 
£>-sphere of radius R. We adopt the convention that a 
D-spheie is the surface of a (D + l)-dimensional ball. 
(Thus, for example, the Berry system is 2-spherium.) 

Quantum mechanical models for which it is possible 
to solve exactly for a finite portion of the energy spec- 
trum are said to be quasi-exactly solvable |81| and we 
have recently discovered that Z?-sphcrium is a member 
of this small but distinguished family [71] [73] . We have 
found that the Schrodinger equation for _D-spherium can 
be solved exactly for a countably infinite set of R values 
and that the resulting wave functions are polynomials in 
the interelectronic distance r 12 = |fi — r 2 |. 

The zeroth-order Hamiltonian of Z?-spherium is 



H 



d9 2 



{D - l)cot6» 



d 

d~6' 



(39) 



where 8 is the interelectronic angle and the associated 
cigenfunctions and eigenvalues are, respectively, 



TABLE II. E (0) , E w , £ (2) and E^l for D-spherium. 



D 



E V>) E W 



E 



(->) 



E 



(2) 



41n2 - 3 
4 368 



3tt 

4 

5 

256 
105vr 

16 

21 
8192 
3465tt 



27tt 2 



64 ln2- 229 
75 375 

24 2 650112 



35 385 875tt 2 

1024, „ 455 803 
In 2 



2205 1 389 150 

4924 588 637011968 



10 395 124 804 708 8757T 2 



Using the partial-wave expansion of r 12 , one finds 

(n + 1) D _2 



12 



where (a) n is the Pochhammer symbol (64) 

T(a + n) 



(a) n 

From this, one can show that 
£7^(D,Sp) = 



T(a) 



(ASp) -r(i.-|)r(f) 



(43) 

(44) 

(45) 
(46) 



and the second-order energy is given by 



E^(D,Sp) = 



n=l 



■ _D-1 N 
, 2 / 



E 



47T 

{n + l) D -2 



(47) 



1^2 
D-l 



D-l 



which reduces to a generalized hypergeometric function. 
It is also easy to show [70] that E^^D, Sp) = 0. 

The E^ 2 ' (and thus E c ) value for 2-spherium was first 
reported by Seidl |79j but elementary expressions for any 
D can be obtained from Eq. ( 47 ) and these are reported 
for 2 < D < 7 in Tableful 



*„(0)=AfC„ 2 (cosfl), 
E n = n(n + D — 1), 

where C n 2 is a Gegenbauer polynomial [64] and 



N = 



l 2 D - 3 {2n + D - l)r {^f 1 ) 2 r(n + 1) 
nT{n + D- 1) 



(40) 
(41) 



(42) 



C. Hookium 

Hookium (Ho) consists of two electrons that repel 
Coulombically but are bound to the origin by the har- 
monic potential 



V(r) 



(48) 
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TABLE III. E {0) , E w , E (2) and E$ for D-hookium. (G is Catalan's constant 



D 


E (0) 


E W 








2 


2 




2G-1T In 2 


-^3(l,4|2,2,2 


••i) 



3 3 



- l--(l + ln2) - - - \l + ln(8 - W3)l 

7T 7T 3 7T L J 



1 4-tt 2G-7rln2 7T „ / , 3 3 „ „ , 

4 4 -J- 1 4F3 1,1, -, -;2,2,3; - 

2V 2 16 4 256 V 2 2 

2/2^ 5 8 „, „ N 88 



3 /tt 104- 27tt 9(2G-7rln2) 3vr ^ / , 3 3 „ „ „ 1 

6 6 8V2_ —512- + ^4 ' -2048 4F3 ( M '2'2 ;2 ' 2 ' 4; 4 

2 89 128 ,„„ „, „ N 416 16 



This system was introduced 50 years ago by Kestner and 
Sinanoglu [82. and solved analytically by Santos [83] and 
Kais et al. |84j for a particular value of the harmonic 
force constant. Later, Taut showed that it is quasi- 
exactly solvable in that its Schrodinger equation can be 
solved for a countably infinite set of force constants [85] . 
An interesting paper by Katriel et al. discusses simi- 
larities and differences between the hookium and helium 
atoms [86] . 

The one-electron Hamiltonian in D-hookium is 



h = 



dr 2 



D — 1 d 

r dr 



2 ' 



and the zeroth-order wave functions are 

D 

#e{ri,r 2 ) = ]J ipa k {xi, k )ip bk (»2,fc), 

k=l 



(49) 



(50) 



where Xi. k is the fcth Cartesian coordinate of electron i, 
and a k and b k are non-negative integers. The orbitals are 
the one-dimensional harmonic oscillator wave functions 

i/> a (x) = \f¥a^T 2 H a {x) exp(-a; 2 /2), (51) 

where H a is the ath Hermite polynomial |64j . The energy 
differences between the eigenstates are given by 



D 

Si- e = ^(a fc + b k ) = 2n, 

k=l 



(52) 



where In is the excitation level, i. e. the number of nodes 
in ^>g. It is not difficult to show that 



£ (0) (£>,Ho) = D, 



E (1) (D,Ko) 



1 r(^i) 



V2 r(f ) 



(53) 
(54) 



Both and J^f 

can be found by direct summation 
|87j . as in Eq. (47). The sum includes all single and 



double excitations for E^ 2 \ but only singles for i%F- The 
integral \ r i2\ ^t) vanishes unless all of the a k + b k 
are even and, in that case, it is given by 



1 r(^i)n 



V2tt T(n + 2) 



D 

n 

k=l 



;a k —b k 



a k + b k + 1 



(55) 



In this way, one eventually finds 

r 



£ (2) (D,Ho) = - 



(ill 

4r(f) 2 fx^)n 



' D-l 
, 2 



2 00 

V 



' D-l 
> 2 



, 2 00 



;d! (1/4)- 



2r (t) "=1 ( 2 'n 



(56) 
(57) 



which reduce to generalized hypergeometric functions. 

.E( 2 )(3,Ho) has been derived by several groups [87l - 
155] . and the energies for other D have been reported in 

Ref. [SS]. Closed- form expressions for E^ and E^,, for 
2 < £> < 7, are listed in Table [TTTI 



D. 



Ballium 



Ballium (Ba) was first studied by Alavi and co-workers 
[90H92] and consists of two electrons, repelling Coulom- 
bically, but confined within a ball of radius R = 1/Z. It 
has been used for the assessment of density-functional 
approximations [5T] H3 [M] and the study of Wigner 
molecules [36] at low densities [92J, [95] [96] . We recently 
obtained near-exact energies for various values of R [9"T] . 

The one-electron Hamiltonian for D-ballium is 



h = - 



dr 2 



D-l d 
r dr 



V(r), 



(58) 
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and the external potential is defined by 
V(r) 



0, if r < R, 
oo, otherwise, 



or equivalently 



V(r) = (r/R) m , m 



— > oo. 



(59) 



(60) 



Any physically acceptable eigenfunction of ^ must sat- 
isfy the Dirichlet boundary condition 



*(ri = R) = *(r 2 = R) = 0. 



(61) 



The associated zeroth-order wave function of the zeroth- 
order Hamiltonian ( 58 ) is 



*o(ri,r 2 ) 



Jd/2(k) 



Jpft-iWl) Jp/2-i{nr 2 ) 

D/2-1 D/2-1 ' V" > 



In (62), k = 3d 12-1,1 and j^,k is the A;-th zero of the 
Bessel function of the first kind J D/2 -i [64 . The 
values are easily obtained from the relation 



£(°)(D,Ba) = k 2 



(63) 



For odd D, E^ can be found in closed form via Eq. (pi 
For example, for D = 3, 



£ (1) (3,Ba) = 2 



1 - 



Si(27r) Si(47r) 



2tt 



47T 



where Si is the sine integral function 
Using the basis functions 



i> n i m = (l-x 2 )(l-y 2 )x 2n y 2l z m , 



where 



x = ri/R, y = r 2 /R, 



z = r 12 /R, 



(64) 



(65) 



(66) 



and n, / and m are non-negative integers, one finds that 
the second-order energy is given by pi]) where 



A = T-£ (0) S, 



b = C J 



£ (1) S - u 



(67) 
(68) 



The vector C contains the coefficients of the zeroth-order 



wave function ( 62 1 expanded in the basis set ( 65 1 . 



The integrals needed to compute the different matrix 
elements are of the form 



T nlm = I .c n y l z m <lr. 
with the volume element 



(69) 



dr = xy z J 2 dxdy dz, (70) 
J = (x + y + z)(x -y + z)(x +y- z)(x -y- z), (71) 



and domain of integration 



1 rl rx+y 
dr = I dx dy dz. 
JO J\sc-y\ 



(72) 



One eventually finds 

p t D—l \ j^n+l+m+2D 



r(f) n + l + m + 2D 



(c+in, (73) 



and 



b _ 3 



I" = 



P ( a+D _t b+D-2 . 
r2 \ 2 ' 2> 2 ' 



D+2 D . 



a + D 



(74) 



(2) 

E^p values can be found using the Byers-Brown- 



Hirschfelder equations (see Sec. II A) and numerical val- 
ues of E^ and E^ are listed in Table [i] 

III. LARGE-DIMENSION LIMIT 
A. The conjecture 

In the large- D limit, the quantum world reduces to 
a simpler semi-classical one |98j and problems that defy 
solution in D = 3 sometimes become exactly solvable. In 
favorable cases, such solutions provide useful insight into 
the D = 3 case and this strategy has been successfully 
applied in many fields of physics [Ml 1100) . 

Following Herschbach and Goodson |101[ 1102) , we ex- 
pand both the exact and HF energies with respect to 
D. Although various asymptotic expansions exist |103j 
for this dimensional expansion, it is convenient [100] to 
write [IU5HTP6] 



E {2 \D.y) 



EW(D,V) 
where 



E (2,o)(y) ^ E^(V) 
D 2 

n(2,0) 



D 2 
E^\V) 
D 2 



L> 3 

2,1 



L> 3 



E?' 1] {V) 



E^Xv) = E^{V)-E (2 f(V) 
E^ 1 \V) = E^ 1 \V)-E < ^\v) 



(75) 
(76) 
(77) 

(78) 
(79) 



Such double expansions of the correlation energy were 
originally introduced for the helium-like ions, and have 
led to accurate estimations of correlation |1071 1108] and 
atomic energies [1091 1110) via interpolation and renor- 
malization techniques. Equation (77) applies equally to 
the 1 S ground state of any two-electron system confined 
by a spherical potential V(r). 

For helium, it is known plU QM1 [m] that 



£^°>(He) 



E^Qie) 



163 

384' 



(80) 
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and we have recently found [S3] that 








1 

~ ~8' 


E^(Ro) = 


111 

ZOO 


(81) 


^ 2 ' 0) (s P ) 


1 

= ~8' 


£f 1} (s P ) = 


53 
~128' 


(82) 


£f°)(Ba) 


1 

~~8' 


^ 2 ' 1) (Ba) = 


85 
~L28' 


(83) 



The fact that Ec 2 ' ^ is invariant to the external potential 
and Ec ' depends only weakly on it explains why the 
high-density correlation energies (Table |l| of all the sys- 
tems are similar, though not identical, for D = 3 [ST)ll97| . 
On this basis, we conjectured [59] that 



EW(D,V) 



8D 2 



con 

L> 3 



(84) 



holds for any spherical confining potential, where the co- 
efficient C(V) varies slowly with V(r). 



B. The proof 

In this Section, we will summarize our proof of the con- 



jecture (84). More details can be found in Ref. [HO]. We 



prove that E^ 2 ' ^ is universal, and that, for large D, the 
high-density correlation energy of the 1 S ground state of 
two electrons is given by ( 84 ) for any confining potential 
of the form 

V(r) = sgn(m)r m v(r), (85) 
where v(r) possesses a Maclaurin series expansion 



v(r) = v + v\r + v 2 — + 



(86) 



After transforming both the dependent and indepen- 
dent variables [60J, the Schrodinger equation can be 
brought to the simple form 



if + U + V + ^W ) <\>l, fjy]' D . (S7) 



in which, for S states, the kinetic, centrifugal, external 
potential and Coulomb operators are, respectively 



-2T = 



o 2 



if 



dr\ dr 2 



U = 



1 



1 



dd 2 



1 



2 sin 2 9 



V 



w 



V(n) + V(r 2 ), 
1 



\A"i + r 2 — 2rir2 cos 9 ' 
and the dimensional perturbation parameter is 
A= (£>-2)(£>-4) _ 



(89) 
(90) 
(91) 

(92) 



(2,0) E (2,0) 



TABLE IV. E (2 ' 0) , 
various systems and v(r) = 1. 



and E, 



(2,1) 



coefficients for 



System 


m 


__g(2,0) 


P (2,0) 


-E {2 ' 0) 


-E^ 


Helium 


-1 


5/8 


1/2 


1/8 


0.424479 


Airium 


1 


7/24 


1/6 


1/8 


0.412767 


Hookium 


2 


1/4 


1/8 


1/8 


0.433594 


Quartium 


4 


5/24 


1/12 


1/8 


0.465028 


Sextium 


6 


3/16 


1/16 


1/8 


0.486771 


Ballium 


CO 


1/8 





1/8 


0.664063 



In this form, double perturbation theory can be used to 
expand the energy in terms of both l/Z and 1/A. 

For D — oo, the kinetic term vanishes and the elec- 
trons settle into a fixed ("Lewis") structure |101| that 
minimizes the effective potential 



X =U + V 



1 



-W. 



The minimization conditions are 

dX(n,r 2 ,9) _ dX( ri ,r 2 ,e) 



dr 2 



= 0, 



dX(r u r 2 ,( 
89 



0, 



(93) 

(94) 
(95) 



(96) 



and the stability condition implies m > —2. Assuming 
that the two electrons are equivalent, the resulting exact 
energy is 

It is easy to show that 

a 2 ( 1 . m + 1 vi \ 1 

' A - + . 

m vq J Z 

a 1 



a - 



+ 2 V2V2 

COS 6*00 = - 



4V2Z 



(97) 
(98) 



where a~( m+2 ^ = sga(m)mvo. 

For the HF treatment, we have = tt/2. Indeed, the 
HF wave function itself is independent of 9, and the only 
9 dependence comes from the Z3-dimensional Jacobian, 
which becomes a Dirac delta function centered at 7r/2 as 
D — > 00. Solving (94 1, one finds that r^ F and are 



equal to second-order in l/Z. Thus, in the large- D limit, 
the HF energy is 



£* F = X 



r HF HF 

' 00 > ' oo > 



(99) 



and correlation effects originate entirely from the fact 
that #oo is slightly greater than 7r/2 for finite Z. 

Expanding ( 96 ) and ( 99 1 in terms of Z and D yields 



EM><y) = -l- 



1 



E 



(2,0) 
HF 



(V) 



8 2(m + 2)' 
1 



2(m + 2)' 



(100) 
(101) 



8 



Coefficients 



-1 £(2,0) 



2 .•" 3 4 5 6 



m 



F (2,0) 
fc HF 



£ (2,0) 




FIG. 2. Coefficients of the exact (dashed), HF (dotted) and 
correlation (solid) energies with respect to m, for v(r) = 1 
(Eqs. f75|, (f76l) and fF7}). 



thus showing that both ij( 2 <°) and depend on the 

leading power to of the external potential but not on v(r). 
Subtracting these energies yields 



E tm {v) 



(102) 



and this completes the proof that, in the high-density 
limit, the leading coefficient E^ 2 ' ^ of the large- £> expan- 
sion of the correlation energy is universal, i. e. it does not 
depend on the external potential V(r). 



The result ( f 02 ) is related to the cusp condition [If 21 - 

eh 



dr 12 



1-12=0 



D - 1 



^D{ri2 = 0), 



(103) 



which arises from the cancellation of the Coulomb oper- 
ator singularity by the /^-dependent angular part of the 
kinetic operator [2]. 

The E^V and E^f 1 coefficients can be found by con- 
sidering the Langmuir vibrations of the electrons around 
their equilibrium positions [lOf 1 If 02] . The general ex- 
pressions depend on vq and V\, but are not reported here. 



However, for v(r) — I, which includes many of the most 
common external potentials, we find 



85 
128 



9/32 



(TO + 2) 3 / 2 

1/2 



1/16 



(m + 2) 1 /2 (to + 2)V2 + 2 ! 



(104) 



(2 1) 

showing that Ec , unlike E, 

(2 1) 

Numerical values of Ec are reported in Table 
various systems, and the components of the correlation 
energy are shown graphically in Fig. [2j 
IV. CONCLUSION 

In this paper, we have reviewed our recent progress in 
the determination of the high-density correlation energy 
for four two-electron systems: the helium atom (He), the 
Hooke's law atom (Ho), two electrons confined on the 
surface of a sphere (Sp), and two electrons trapped in a 
ball (Ba). In the large- Z limit for D = 3, we have found 



n(2,0) 



, is potential-dep ende nt. 

for 



IV 



£ c (He) ps £ c (Sp) ps £ c (Ho) « £ c (Ba). 



(105) 



These striking similarities can be rationalized by treating 
the dimensionality D of space as a system parameter, and 
we have proved that, as D grows, all such correlation 
energies exhibit the same universal behaviour 



-1/(8C 2 



(106) 



in a D-dimensional space. This is true irrespective of the 
nature of the external potential that confines the elec- 
trons. 
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